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ABSTRACT 

Context. Within the core accretion scenario of planetary formation, most simulations performed so far always assume the accreting 
envelope to have a solar composition. From the study of meteorite showers on Earth and numerical simulations, we know that 
planetesimals must undergo thermal ablation and disruption when crossing a protoplanetary envelope. Thus, once the protoplanet has 
acquired an atmosphere, not all planetesimals reach the core intact, i.e., the primordial envelope (mainly H and He) gets enriched in 
volatiles and silicates from the planetesimals. This change of envelope composition during the formation can have a significant effect 
in the final atmospheric composition and on the formation timescale of giant planets. 

Aims. To investigate the physical implications of considering the envelope enrichment of protoplanets due to the disruption of icy 
planetesimals during their way to the core. Particular focus is placed on the effect on the critical core mass for envelopes where 
condensation of water can occur. 

Methods. Internal structure models are numerically solved with the implementation of updated opacities for all ranges of metallicities 
and the software CEA to compute the equation of state. CEA package computes the chemical equilibrium for an arbitrary mixture of 
gases and allows the condensation of some species, including water. This means that the latent heat of phase transitions is consistently 
incorporated in the total energy budget. 

Results. The critical core mass is found to decrease significantly when an enriched envelope composition is considered in the internal 
structure equations. A particular strong reduction of the critical core mass is obtained for planets whose envelope metallicity is larger 
than Z a 0.45 when the outer boundary conditions are suitable for condensation of water to occur in the top layers of the atmosphere. 
We show that this effect is qualitatively preserved also when the atmosphere is out of chemical equilibrium. 

Conclusions. Our results indicate that the effect of water condensation in the envelope of protoplanets can severely affect the critical 
core mass, and should be considered in future studies. 

Key words. Planet formation, exoplanet atmospheric composition. 


1. Introduction 


The metallicity of the giant planets of our solar system shows 
that they are considerably enriched in heavy elements compared 
to solar composition (Guillot & Gautier 2014, and references 
therein). Even though this change of composition with respect to 
solar could be a consequence of an eroding core (Guillot et al. 
2004; Wilson & Militzer 2012), numerical studies have shown 
that the change of composition induced by the disruption of plan¬ 
etesimals in the envelope of growing giant planets can be consid¬ 
erable (Podolak et al. 1988; Mordasini et al. 2006; Iaroslavitz & 
Podolak 2007; Fortney et al. 2013). For instance, a protoplanet 
with a core of 6 M@ can have an envelope dense enough to com¬ 
pletely destroy icy planetesimals of up to 100 km before they 
reach the core (Baraffe et al. 2006). The dominant size of plan¬ 
etesimals at the time planet formation takes place might even 
be smaller, as is suggested by studies which incorporate the col- 
lisional fragmentation in the evolution of planetesimals (Inaba 
et al. 2003; Kobayashi et al. 2011; Guilera et al. 2014). This 
means that a considerable enrichment can readily occur even for 
protoplanets with small cores. Thereby, the contribution of met- 
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als (compounds heavier than H and He) to the total mass of the 
envelope of a protogiant planet starts to be relevant already in 
the early stages of growth. We show in this paper that taking into 
account this early enrichment of the envelope can have signifi¬ 
cant consequence on the subsequent growth and evolution of a 
planet. 

The aforementioned studies (Podolak et al. 1988; Mordasini 
et al. 2006; Iaroslavitz & Podolak 2007) have focused on com¬ 
puting the trajectory of planetesimals during their way to the 
core, as well as their mass and energy deposition in the envelope 
due to thermal ablation and disruption. However, the resultant 
metallicity was not included in the thermochemical properties 
of the envelope, i.e., it was not taken into account to solve the 
internal structure equations. 

The opposite focus was given by Hori & Ikoma (2011) (here¬ 
after HI11). In the first part of their work, they assumed a given 
envelope metallicity and incorporated a consistent equation of 
state and opacity tables to solve the static internal structure equa¬ 
tions. “Static” means that the only source of energy considered is 
the release of gravitational energy due to the accretion of plan¬ 
etesimals. This implies that the energy release due to the con¬ 
traction of the envelope (which necessarily provides a varying 
luminosity with time) is neglected. 
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Static internal structure calculations allow for the computa¬ 
tion of the critical core mass, which corresponds to the max¬ 
imum mass of the core able to sustain an envelope in hy¬ 
drostatic equilibrium (Perri & Cameron 1974; Mizuno et al. 
1978; Mizuno 1980). Analytical work by Stevenson (1982) and 
Wuchterl (1993) (for fully radiative and fully convective en¬ 
velopes, respectively) show that when the core becomes criti¬ 
cal, the mass of the envelope is of the same order of magnitude 
than the core mass. This means that at the critical point the self¬ 
gravity of the envelope becomes important. Further, quasi-static 
studies of planet growth (Bodenheimer & Pollack 1986; Pollack 
et al. 1996), show that when the energy release due to the con¬ 
traction of the envelope is taken into account, rapid gas accretion 
takes place when M core ~ M env , that is, when the core is around 
criticality. Since the timescale of gas accretion at critical mass 
is much shorter than the one required to build a critical core, 
reaching criticality is a synonym of being able to form a giant 
planet. Hence, the relevance of the critical core mass concept in 
the core-accretion scenario. 

The enrichment of the envelope due to planetesimal disrup¬ 
tion causes a reduction of the critical core mass due to the in¬ 
crease in mean molecular weight of the envelope (first pointed 
out by Stevenson 1982), but more importantly, due to the reduc¬ 
tion of the adiabatic temperature gradient that takes place when 
chemical reactions among volatile molecules occur (Hill). In 
this work we have included the effect of water latent heat and 
studied its effect on the critical core mass. 

The role played by the condensation of water (i.e., the for¬ 
mation of clouds) in Earth’s atmosphere is well known among 
the climate community (see, e.g., Pierrehumbert 2010, chap. 2). 
When vapor raises above the condensation level, water droplets 
form and in the phase transition process latent heat is released. 
This modifies the energy budget of the atmosphere, causing the 
pressure-temperature profile to switch from dry to moist adiabat. 
We will show that this change in the temperature profile actually 
can lead to a severe reduction of the critical core mass. 

The paper is organized as follows. The basic equations, as¬ 
sumptions and numerical scheme are described in Sect. 2. The 
opacities and equation of state implemented are described in 
Sect. 3 and 4, respectively. Sect. 5 is devoted to the results con¬ 
cerning the effect on the critical core mass due to water conden¬ 
sation in an envelope of uniform metallicity. In Sect. 6 we justify 
our assumptions concerning the use of an ideal equation of state 
and uniform metallicity. We discuss as well the dependence of 
our results on the boundary conditions, and the implications of 
our results on the diversity of planets that can be formed. Finally, 
Sect. 7 summarises our main conclusions. 


2. Numerical approach 

We have developed a code that solves the static internal struc¬ 
ture equations assuming spherical symmetry and a constant lu¬ 
minosity (or constant accretion rate of solids when specified). 
For better numerical performance, the equations are written as 
a function of v = r 3 /M r , where r is the radial coordinate of the 
planet, and M r the mass enclosed inside a sphere of radius r. 
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where P is the pressure, p, ol the total density and p gas the gas 
density , T the temperature, and G the gravitational constant. 
The gradient in Eq.(3) is defined as V s For the radiative 

case it is equal to: 

_ 3 kLP 

rad " 64no-GMJ A 


cr being the Stefan-Boltzmann constant, k the opacity and L the 
luminosity. For the convective transport V is the gradient evalu¬ 
ated at constant entropy (adiabatic gradient): 



Convection occurs if V a[ j < V rac j (Schwarzschild criterion). 


The three differential equations are solved using a 4th- order 
Runge-Kutta algorithm. Two boundary conditions are specified 
at the top of the gaseous envelope (7', )U , and P out ) and one at the 
core (core mass or core density). The outer boundary conditions 
and the position of the planet are listed in Table 1. The total mass 
of the planet (Mp) is defined as the mass inside the Bondi or Hill 
radius, whichever is smaller. They are, respectively: 



R h = 


I Mp \ 

13 mJ 


1/3 


( 6 ) 


where c s is the sound speed and a the semimajor axis of the 
planet. 

Two different schemes to solve the above equations were im¬ 
plemented. In the first approach (used for the results shown in 
Fig.8), the core mass (M core ), the core density (p c ore ), and the 
accretion rate of solids (M c ore) are given. Then the total mass 
of the planet (Mp) is computed through an iterative process. A 
first value for the total mass is assumed and the structure equa¬ 
tions are solved starting from the outer boundary and stopping at 
the core. The total mass is iterated until the mass of the core is 
reached within a given tolerance. Since in this approach the core 
radius (R core ) and core mass (M core ) are known a priori, the lumi¬ 
nosity of the envelope is given by the release of the gravitational 
potential energy of the planetesimals: 

T _ GM core M core 4 

L ~ R 

A. core 

which implies that the energy is assumed to be deposited at the 
surface of the core. This is not consistent with our assumption of 
the disruption of the planetesimals occurring before they reach 
the core, but to do this properly, a varying luminosity in the en¬ 
velope should be implemented, computed from the mass and en¬ 
ergy deposition of the disrupted planetesimals. This will be done 
in a future work. However, for the cases of highly enriched en¬ 
velopes (the focus of this work), the envelopes are fully convec¬ 
tive (see Sect.6.1), and hence, the thermal structure is indepen¬ 
dent on L. Thus, the precise value of L is not expected to affect 
our results. 

In the second approach, the total mass of the planet and total 
luminosity are given. M core is not known a priori because M core is 
not known. In this case the integration is also performed from the 
surface to the core, but the difference is that now Mp is known, 

1 p to , includes the mass of all species and p gas just the one of the 
gaseous component. These two densities are equal if no condensation 
takes place. See Sect. 4.2 for more details. 
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a 

5.2 AU 

Tout 

150 K 

Tout 

0.267 dyn/cm 2 

P core 

3.2 g/cm 3 

L 

1.0 x 10 27 erg/s 


Table 1: Parameters used in all simulations. The planet location 
and boundary conditions defined here are referred in the main 
text as “standard case”. These parameters are only modified in 
Sect. 5.4. 


so no iteration is needed. The integration stops when p r (p r = 
is equal to the density of the core. The mass inside this 
radius is then defined as M core . The independent variable v in 
Eqs. 1,2, 3 was introduced with the purpose of making it possible 
for the integration to stop exactly at p r = p core . Otherwise, an 
interpolation in r would be required to obtain the mass at which 
Pr = Pcore ■ This last method is the one used in all the results 
shown in this work, except for those shown in Figs. 2, 3 and 
8. Its advantage resides in the fact that no iterations are needed 
to obtain M c0le , which makes the scheme more stable for core 
masses close to the critical core mass. 


3. Opacities 

Most of the opacity tables available in the literature are for solar 
metallicity. In order to study the effect of enrichment of the en¬ 
velopes of protoplanets in a consistent way, we must use opacity 
tables for all ranges of metallicities. In this work we implement, 
as in Hill, gas opacity tables for all metallicities computed by 
J. Ferguson (based on the calculations of Alexander & Fergu¬ 
son 1994) and Rosseland mean opacities of dust calculated by 
Semenov et al. (2003). 

The Ferguson opacities assume scaled solar abundances (the 
relative abundances of the various elements is assumed solar but 
the total abundance is scaled). The range of validity of the tables 
is for T > 1000 K and -2 < log R < 6, where R — p/T 2 3 , p 
being the gas density and 7/, the temperature in million kelvin. 
Whenever the envelope values of temperature and density lay 
outside these ranges, the tables are extrapolated. Using a linear 
or constant extrapolation does not affect our results. 

As in Hill, we define the total opacity as: 

K — Kgas + ,/Alusl (8) 

where /c gas and /e duM are the Rosseland mean opacity of gas and 
dust grains, respectively. The reduction/enhancement factor / in¬ 
cludes the unknown values of the size of the dust grains in proto¬ 
planetary envelopes. In principle, the grains released by the abla¬ 
tion of planetesimals could increase the dust opacity with respect 
to the disc values, but also the grains could rapidly coagulate and 
hence increase their size, reducing in this way the opacity of the 
dust (Podolak 2003; Movshovitz & Podolak 2008; Movshovitz 
et al. 2010; Mordasini 2014; Ormel 2014). 

In this work we adopt values of / = 1, since the effect of 
other reduction/enhancement values was already investigated by 
Hill. Using / = 10 hardly affects the critical core values with 
respect of using f — 1, and using lower values of / reduces the 
critical core mass even further. Thereby, setting / = 1 provides 
upper values for the critical core mass as far as the effect of dust 
opacity is concerned. Regardless of the fact that these are upper 
limits, we will show that they become increasingly smaller with 
increasing metallicity. 


In Sect.6.2.2 we make use of the Freedman et al. (2008) 
opacities, which are detailed computed opacity tables for H and 
He and solar metallicity. We use these opacities to illustrate the 
effect of different opacity tables in the case of envelopes with ex¬ 
tended radiative regions. All other results were computed mak¬ 
ing use of the opacity given by Eq.8 . 

4. Equation of state 

Beyond suitable opacities, an enriched envelope needs also to be 
described with a proper equation of state. It turns out that, to the 
best of our knowledge, there is no equation of state available for 
an arbitrary metallicity which takes also into account the degen¬ 
eracy pressure of free electrons. The latter effect play an impor¬ 
tant role in the more massive objects. One EOS accounting for 
these effects is the Saumon et al. (1995) EOS for a mixture of H 
and He (hereafter SCVH). 

In this work we make use of the publicly available software 
CEA (Chemical Equilibrium with Applications) developed by 
NASA (Gordon et al. 1994). This program solves the chemical 
equilibrium equations for an arbitrary defined gas mixture. It as¬ 
sumes an ideal gas equation of state for the mixture, but takes 
into account the proper dissociation and ionization temperatures 
of each compound formed. Despite the fact that CEA does not 
consider the pressure of degenerate electrons (which translates 
into an overestimate of central temperature and pressure for mas¬ 
sive objects), we have used it to study the effect of condensation 
on the envelope structure. 

CEA has the advantage of containing a database of ~ 1000 
compounds, and the fact that it considers condensed species, like 
HiO in solid and liquid phase. 2 The latent heat of phase transi¬ 
tions is included in the calculation of the thermodynamical quan¬ 
tities. For an explanation on the calculation of the density when 
condensation takes place, see Sect.4.2. For a discussion on the 
condensation of other species than water, see Sect. 6.4. 

The program can be customized to use as an input the pres¬ 
sure, temperature and composition in molar or mass fraction. In 
our simulations, the envelope is assumed to be made of hydro¬ 
gen, helium, carbon and oxygen. We take a fixed C/O mass ratio 
of 0.24:0.69, which is the one corresponding to Comet Halley 
(Mumma et al. 1993), value that we assume representative of the 
early solar system. The metallicity, Z, is defined as: 


where M/ t en v is the mass of elements heavier than H and He in 
the envelope, and M env the total mass of the envelope. 

Carbon and oxygen constitute the Z-component of the enve¬ 
lope. Z is a free parameter in our study, and is assumed to be 
uniform throughout the envelope. This assumption will be jus¬ 
tified in Sect. 6. The mass abundance of H and He (X and Y, 
respectively) are computed as a weighted average between the 
solar and Comet Halley values, as in HI11. 3 

CEA gives as an output the density, mean molecular weight, 
internal energy, specific heat, and logarithmic derivatives of cer¬ 
tain thermodynamical quantities from which it is possible to 
compute the adiabatic gradient. 

2 For the complete list of condensed species consid¬ 
ered in CEA, see Appendix B of the User’s Manual from 
http://www.grc.nasa.gov/WWW/CEAWeb/ceaWhat.htm 

3 This choice of planetesimal composition is dictated by the compar¬ 
ison with earlier work of HI11. These ratios could be varied to match, 
e.g, the composition of different stars. 
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The adiabatic gradient (defined by eq. 5) is an important 
quantity because, together with the radiative gradient (eq. 4), it 
determines the type of heat transport present in the envelope. 
Given its relevance to our study of the structure of planetary en¬ 
velopes, we present in the next section how we determine this 
gradient. 

4.1. Adiabatic gradient with CEA 

One of the outputs of CEA is the log-derivative of the volume 
with respect to temperature at constant pressure. From this quan¬ 
tity and the specific heat at constant pressure, the adiabatic gra¬ 
dient can be easily calculated. 

From the implicit function theorem: 



And the gas density must satisfy: 4 

_ MgasP 
Pgas ~ rj. 

The gas density is readily obtained from: 


(18) 


Pgas — (1 Xc)piot (19) 

This density enters in the hydrostatic equilibrium equation, be¬ 
cause the pressure gradient is exerted just by the gas . 5 For the 
mass conservation equation, the density that includes the total 
mass, i.e. p tol , is the one that must be used. (See Sect.2). 


5. Critical core mass for different uniform 
metallicities 

5.1. Comparison with HI11 


From the definition of the adiabatic gradient, 




Then, 


P(dS/dP) T 

T(dS/dT) P 


( 11 ) 


( 12 ) 


From the definition of the heat capacity at constant pressure 

(Cry. 



(13) 


Deriving the Gibbs free energy (G - U - 'I S + PV) with respect 
to T and P it can be shown that: 



(14) 


And finally, substituting by the log-derivatives and using the fact 
that the mixture is an ideal gas: 


^ nR (dlnV\ 
ad = — p \iMr)p 


(15) 


where n is the number of gas moles and R the gas constant. 


4.2. Condensed species in CEA: density calculation 

In CEA the mean molecular weight is defined as M — — (M to , 

n gas 

being the total mass of the system and n gas the number of moles 
of the gas), which is not the mean molecular weight of the gas in 
the case when condensation occurs. 

It can be shown that the molecular weight of the gas satisfies: 


Mgas = (1 - X C )M 


(16) 


where X c is the mass fraction of the condensed species, quantity 
readily computed with CEA using Eqs. 2.3b and 2.4a of Gor¬ 
don et al. (1994). Since the density computed by the program 
satisfies: 


_ MP 

Ptot — 


(17) 


The static internal structure equations were solved for a range 
of uniform envelope metallicities between 0.1 and 0.95, and a 
wide range in protoplanetary mass. The planet location, bound¬ 
ary conditions and envelope composition were taken as in HI 11 
(see Table 1), in order to be able to perform comparison tests. 
The equality in composition means that we set the same X and 
Y for a given Z (see Sect.4), and that we allowed the same 13 
species in gaseous phase (H, He, O, C, H2, O2, CO2, H2O, CH 4 , 
H + , O and e ) to be formed when chemical equilibrium is com¬ 
puted. Hereafter we will refer to this case as the “restricted” case. 

A comparison with HI11 is shown in Fig. 1 (green versus red 
curve). The agreement between the two calculations is remark¬ 
able, especially considering that a different way of computing the 
EOS was used in the different works (HI11 did not use CEA, but 
computed the thermodynamical quantities partly from the NIST 
tables, see HI11). 

5.2. Inclusion of H 2 0 condensation 

The blue line of Fig. 1 shows the results for the same setup, but 
in this case no restriction is imposed in the compounds that can 
be formed when the chemical equilibrium is computed. Hence, 
CEA makes use of its chemical database, and allows ~ 50 com¬ 
pounds to be formed from the same proportion of H, He, C and 
O than before. It also allows condensed species to be formed. 
We will refer to this case, where no restriction is imposed in the 
species that can be formed at equilibrium, the “non-restricted” 
case. 

Regardless of the case, the most remarkable feature shown 
by Fig.l is the drastical reduction of the critical core mass with 
increasing Z. Comparing now the “non-restricted” with the “re¬ 
stricted” case, it can be seen that the reduction of M cr ; t with in¬ 
creasing Z is larger for the first case. Namely, for Z = 0.7, the 
difference in critical core mass between the restricted and the 
non-restricted case, is of AM cr i t ~ 2M 0 , and since the critical 
mass of the non-restricted case is M cr ; t ,NR ~ 0.2M@, the relative 
difference is t Met " « 10 . 

M cr it,NR 

4 or equivalently, p gas = being m A the atomic mass unit and k B 

the Boltzmann constant. The mean molecular weight of the gas, p gas , is 
numerically equivalent to M gas but is unitless. 

5 It is debated in the literature if actually p gas or p, ot should be used 
in the hydrostatic equilibrium equation, but recent work suggest that 
the correct expression is with p gas (Gorshkov et al. 2012). In any case, 
we have checked that using p gas or p tot in the hydrostatic equilibrium 
equation does not affect our results. 
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Fig. 1: Critical core mass as a function of metallicity. Green line: 
results of the case “wholly polluted” of Fig.2 of Hill, red line: 
this work, restricted case. Blue line: this work, non-restricted 
case. (See main text for the definitions of the different cases). 


In order to understand the difference in M cm between the 
“restricted” and “non-restricted” cases, we plot the gas mean 
molecular weight and adiabatic gradient as a function of the en¬ 
velope temperature for Z = 0.7. For both profiles we assume the 
same core mass, which corresponds to the critical core of the 
non-restricted case for Z = 0.7. The purpose of this choice is 
to ensure the same core gravity. This, in turn, ensures that the 
differences observed in the atmospheric profiles originate from 
the chemistry and not from a difference in initial pressure gra¬ 
dients. The results are shown in Figs. 2 and 3. A decrease of 
the gas mean molecular weight and of the adiabatic gradient in 
the outer layers are observed when comparing the non-restricted 
with the restricted case (solid and dotted lines, respectively). By 
analyzing the composition in the outer layers of the envelope of 
the non-restricted case, we find that the 13 species assumed by 
HI11 are indeed the dominating ones in the CEA calculations as 
well, as shows Fig. 4 for Z = 0.7. 6 However, we also find that 
large amounts of water are being condensed in the coldest layers 
of the envelope of the non-restricted case, fact which does not 
occur in the restricted case (see Fig. 4). In fact, observing Figs. 
2 and 4, it can be appreciated that the remarkable difference in 
mean molecular weight in the outer part of the envelope between 
the non-restricted and restricted cases lasts as long as H 2 0 coex¬ 
ist in both vapor and solid phase. This is a consequence of what 
was stated in Sect.4.2: the larger the mass fraction of condensed 
species, the lower p gas (see Eq.16). 

Water condensation in the outer layers of the envelope mod¬ 
ifies the molecular weight and adiabatic gradient. In order to as¬ 
sess which of the two effects plays a role in reducing M CI i t , we 
performed the following test. We computed the envelope profile 
for Z = 0.7 and Mp = 0.29 M e using p of the restricted case and 
V a d of the non-restricted case. The profile follows exactly the 
same profile as the non-restricted case. Therefore, the reduction 
of the adiabatic gradient is what causes the smaller critical mass. 

A straightforward way to understand why the reduction of 
the adiabatic gradient tends to lower the critical core mass is as 
follows: given a fixed core mass (as in the cases shown in Fig.3), 

6 The envelope composition for Z = 0.2 and 0.4 are also shown in 
Fig. 4, to illustrate how the abundance of volatiles, in particular water, 
increases with increasing metallicity. 


Fig. 2: Profile of the gas mean molecular weight of the enve¬ 
lope as a function of its temperature for Z = 0.7 when M core = 
M cr ;t,NR (Z = 0.7). Solid lines: non-restricted case. Dotted lines: 
restricted case. Note that in the the non-restricted case, p gas de¬ 
creases in the outer layers (i.e., low temperature values) when 
comparing to the restricted case. This effect is explained in the 
main text. The difference between both profiles at high tempera¬ 
tures is a consequence of a difference in the pressure-temperature 
profiles, see Fig. 5. 



T [K] 


Fig. 3: Profile of the adiabatic gradient of the envelope as a func¬ 
tion of its temperature for Z = 0.7. The core mass and labels are 
set just as in Fig. 2. Note that there is a strong reduction of the 
adiabatic gradient in the outer layers of the envelope (i.e., low 
temperature values) in the non-restricted case. This effect is ex¬ 
plained in the main text. The difference between both profiles at 
high temperatures has the same explanation as in the previous 
figure. 


the gravity is fixed, and therefore the pressure gradient (Eq. 1) 
is fixed. The thermal profile is established by Eq.(3). From the 
ideal gas equation we have: 


pkT 

pm A 


( 20 ) 


(where p is the mean molecular weight of the gas mixture and m A 
the atomic mass unit). If the envelope is convective and V a d di- 
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minishes, the temperature for a given radius diminishes. Hence, 
in order to keep the same pressure gradient, the density of the 
envelope has to increase. Since the envelope mass is the inte¬ 
gral of the density, if the density increases, the envelope mass 
increases as well. The more massive an envelope gets for a given 
core mass, the more likely it is that the core becomes critical (see 
Fig.7). 

As stated before, the formation of other compounds in the 
non-restricted case does not play a role, since the 13 compounds 
assumed by Hill are the most abundant ones when chemical 
equilibrium is computed. But taking into account the conden¬ 
sation really changes the picture, as illustrated in Fig.5. In this 
figure, the pressure-temperature profile for Z=0.7 for the non- 
restricted case, and the restricted case without condensation of 
water, are shown. Superposed, we plot the case where the com¬ 
position is “restricted’', but where water is allowed to condense 
(hereafter we will refer to this case as “restricted-with-water- 
condensation”). The vapor pressure for this last case is also plot¬ 
ted. From this figure it is clear that H 2 O is being condensed from 
the top of the atmosphere until '[’=282 K, since in this range of 
temperature the vapor pressure is equal to the saturation pressure 
of water (green- dotted line, given by the Clausius-Clapeyron re¬ 
lation). As the temperature increases from the top of the envelope 
to the surface of the core, there is a location from which water 
can only exist in vapor phase. For the mass abundance of water 
when Z = 0.7, this happens at T=282 K. If the abundance of wa¬ 
ter is lower than in the Z = 0.7 case, the shift between “moist” 
and “dry” adiabat (as is usually called in atmospheric studies, 
see, e.g., Pierrehumbert 2010, chap. 2) occurs for lower temper¬ 
atures (see Fig.6). This is because the vapor pressure is propor¬ 
tional to the mole fraction of vapor (law of partial pressures). 
This depends on the amount of water present in the outer lay¬ 
ers of the envelope, which increases with metallicity, as shows 
Fig.4. 

Returning to Fig.5, it can be noticed that due to water con¬ 
densation, < (jp)dn- i* 1 l h e coldest layers of the atmo¬ 

sphere. This is the cause of the strong reduction in the adiabatic 
gradient of the non-restricted case shown in Fig.3. 

Fig.5 can also be used to understand why condensation of 
water reduces M cr j t . If we look at the envelope profile from the 
surface of the core to the top of the atmosphere, when we com¬ 
pare the non-restricted with the restricted case, we see that at the 
pressure (equivalent to atmospheric height) where water starts 
to condense, the temperature of the non-restricted case does not 
drop so much as in the restricted case. Since latent heat is re¬ 
leased in the transition vapor-solid and vapor-liquid, the enve¬ 
lope gets warmer when condensation of water takes place. The 
release of latent heat is what really makes the difference be¬ 
tween the “moist” and “dry” profiles. It modifies the pressure- 
temperature profile remarkably, making it possible to have a 
much higher pressure for a given temperature, which translates 
in a much higher envelope density, and therefore, in a consider¬ 
able reduction of critical core mass. 

5.3. Increasing envelope metallicity 

In the previous section, we confirmed an earlier finding of Hori 
& Ikoma (2011) that an enrichment of the envelope leads to a 
reduction of the mass of the core at which gas can be accreted in 
a runaway fashion. In addition, we pointed out that in case the 
enrichment leads to condensation, the reduction of the critical 
core mass is even greater (see Figs.l and 7) and can actually 
become the dominant effect. 



Fig. 5: Pressure as a function of temperature for Z-0.1 and 
Mp = 0.29 M®. In the case where water is allowed to con¬ 
dense (orange lines), the vapor pressure follows, until T=282 
K, the Clausius-Clapeyron curve, which defines the loci where 
the vapor and liquid phase, and vapor and ice phase, can co-exist 
(green-dotted line). Hence, in these range of temperatures, water 
is being condensed, and the pressure-temperature profiles follow 
the “moist adiabat". 



Fig. 6: Pressure as a function of temperature for Z=0.2, 0.4 and 
0.7. Mp = 0.29 M® for all cases. The dotted-black lines corre¬ 
spond to regions where the envelope is radiative. This region is 
extremely tiny in the case of Z=0.2 profile. 


To pose this enrichment problem in an evolutionary context, 
it is interesting to see what would be the total mass of a planet 
that becomes critical from a gradual enriching process, and how 
enriched the envelope becomes when reaching the critical core 
mass. In order to tackle this, we performed the following calcula¬ 
tion. We assumed that a core grows up to 6 M® due to the accre¬ 
tion of planetesimals that reach the core without being destroyed 
in the envelope. In other words, all the planetesimals reach the 
core without undergoing any physical disruption or thermal abla¬ 
tion until the core acquires a mass of 6 M®. Once M core = 6 M®, 
we assume that the protoplanet has acquired an atmosphere suffi¬ 
ciently massive to completely disrupt the infalling planetesimals 
before they reach the core. This choice of 6 M® as a threshold is 
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completely arbitrary. It corresponds to the case of icy planetesi- 
mals with a radius of 100 km or less (Baraffe et al. 2006). How¬ 
ever, as we explained in Sect. 1, planetesimals could be smaller. 
In that case, full disruption could already occur for planets with 
a smaller core; a case also analysed at the end of this section. 

For the case of a fixed core mass of 6 M® as the core thresh¬ 
old for completely disrupting planetesimals in the envelope, 7 
we solved the internal structure equations, and let the envelope 
metallicity grow gradually until the planet becomes critical. The 
results are shown in Fig. 8, where the total mass of the envelope 
as a function of the total mass of the heavy elements acquired 
by the planet is plotted. If initially, when the fully disruption of 
the planetesimals starts, Z = 0.02 (Z G ), then the envelope is able 
to accrete planetesimals until reaching a metallicity of ~ 0.48, 
which represents an envelope mass content of heavy elements 
of 1.4 M 0 . At this point the core becomes critical, i.e., no more 
solution to the static internal equations can be found. Accord¬ 
ing to the standard view of the core-accretion model, the sub¬ 
sequent evolution of the protoplanet would imply the accretion 
of large amounts of gas in a runaway fashion (Mizuno 1980; 
Bodenheimer & Pollack 1986), making it possible for the pro¬ 
toplanet to become a giant planet soon after the protoplanet be¬ 
comes critical. The results found here suggest that after the dis¬ 
ruption of planetesimals in the envelope begins, the planetary 
envelope does not need to grow much in order for the planet to 
become critical. In the particular case of the simulations shown 
here, a protoplanet with M core = 6 M Q needs to accrete 1.4 M® 
of H and He and the same amount of heavy elements to become 
critical. A consistent evolutionary study, where the metallicity is 
computed from the accretion rate of solids, will be done in the 
future to be able to account for the timescales to trigger runaway 
of gas when the envelope enrichment is taken into account. 

A similar analysis of the enrichment process can also be 
drawn for cores smaller than 6 M® by just making use of Fig. 1. 
If we assume, for instance, that the core stops to grow at 1 
M e , then the envelope will get gradually enriched until the core 
becomes critical at Z « 0.62 (see blue curve of Fig.l). The 
corresponding mass of the planet in this critical case is Mp = 
1.5 M®, which means that in this case the mass of heavy el¬ 
ements in the envelope when M cr ; t is reached would be ~ 0.3 
M®, and the mass of H-He ~ 0.2. M®. It would be interesting 
also for cases like this, where ;V/ cn i is low and the envelopes are 
heavily polluted, to study the subsequent evolution. The fate of 
these small but highly enriched critical protoplanets is not easy 
to assess a priori, since the gravity of the core could not be 
large enough to trigger gas runaway (Ikoma et al. 2000; Hori 
& Ikoma 2010). In principle, if runaway of gas accretion cannot 
be reached, the planet would continue its growth by slowly con¬ 
tracting (and therefore, also slowly) accreting gas. If runaway of 
gas is not triggered during the lifetime of the gas disc, the re¬ 
sult could be the formation of a super-Earth, heavily enriched 
planet, which would be a very interesting scenario, given the in¬ 
creasing amount of Earth and super-Earth sized planets that are 
being found (Petigura et al. 2013), and the inference of a very 
enriched composition for these type of objects (Valencia et al. 
2013; Gillon et al. 2012; Nettelmann et al. 2011). A more ex¬ 
tended discussion concerning the different type of planets that 
could be formed when the enrichment of the envelope is consid¬ 
ered, is presented in Sect.6.5. 


7 Since we need to fix the core mass at the beginning of the simu¬ 
lations, we use, in this case, the first numerical scheme described in 
Sect.2. 
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Fig. 7: Total mass of the planet as a function of the mass of 
the core for different uniform metallicities and L = 10 27 erg/s. 
Here the composition of the envelope corresponds to the “non- 
restricted case” (see main text). The highest Z shown is 0.5 for 
aesthetic reasons. 
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Fig. 8: Envelope mass as a function of the mass of solids con¬ 
tained in the planet (Mz). The mass of the core is fixed at 6 M®. 
An accretion rate of solids of 3 x 10~ 6 M®/yr was assumed, which 
correspond to a luminosity of 10 27 erg/s (the same used in all the 
other results). The decrease of M env for low metallicities (shown 
in the insert) is due to the increase in the gas opacity and the 
fact that for these cases of low Z, a radiative zone develops in 
the outer layers. For Z > 0.1, the structures are almost fully con¬ 
vective. Therefore the increase in /c gas does not play a role, the 
dominant effect is the increase in M env due to the inclusion of Z 
in the EOS. 


5.4. Dependence on the boundary conditions and semi major 
axis 

All the results shown before were for a planet whose boundary 
conditions are those stated in Table 1. In order to check the range 
in boundary conditions for which condensation can lead to a sig¬ 
nificant reduction of the critical mass, we perform the following 
test. We set the outer temperature and pressure of the planet for 
a given position in the protoplanetary disc making use of the fol¬ 
lowing simple disc model (see, namely, Armitage 2010): 


Article number, page 7 of 16 















A&A proofs: manuscript no. paper_venturini 



where ciq, Tq, Pq are the standard boundary conditions used in 
all previous results, i.e, those given in Tablel. In this section, 
we will refer to the simulations that use boundary conditions at 
a - ao as the “reference case”. 

The results for the critical core mass as a function of metal- 
licity for different semimajor axis (with the corresponding outer 
temperature and pressure given by Eqs.21, 22) are shown in 
Fig.9. Results for 2 AU < a < 8.5 AU are shown in solid lines, 
superposed with the reference cases in dashed lines (see cap¬ 
tion). The qualitative similarity of the different curves with the 
reference case where water condenses (blue curve) or with the 
one where water is not allowed to condense (yellow curve) pro¬ 
vides a hint on whether water condensation is taking place or 
not for the different semimajor axis considered. We note that 
for a- 2 AU, P out is always too high for water condensation 
to take place. For all the other curves condensation of water 
takes place at some point. The start of the condensation corre¬ 
sponds to the change in the slope of the curves. A particular case 
takes place for a = 3.0 AU, where a valley in M cr i t occurs for 
0.865 < Z < 0.905 (see Fig.9). The drop in A7 cm corresponds 
to the fact that condensation starts to happen at Z = 0.865. The 
same happens for a = 3.2 AU at Z = 0.75 and for a = 3.5 AU 
at Z = 0.55. The surprising thing about the 3 AU case is that for 
even higher metallicities, water condensation ceases. The reason 
for this is as follows: when the minimum in the “valley” of M cl ; t 
is reached, the metallicity has increased so much that there is 
no more free kU at the top of the envelope. That is, IU starts to 
be the limiting reactant, so increasing Z (that is, adding more C 
and O) does not lead to the formation of more water. Instead, 
more CCT starts to be formed in this region of the valley (where 
dM cl -it/dZ > 0). At Z = 0.905, the amount of water formed is low 
enough to make the vapor pressure be lower than the saturation 
vapor pressure for all the P,T values of the envelope. Thereby, 
for Z > 0.905, the profiles follow again the “dry” adiabat, no 
water condensation takes place anymore. 

Despite this particular case of a — 3 AU, the general trend ob¬ 
served in Fig.9 is that lUO condensation occurs more readily for 
planets located further out in the disc. However, while the trend 
is clear, the exact behavior is more difficult to predict. For in¬ 
stance, the behavior of M cr i t changes considerably for 3 < a < 4 
AU, but remains very similar for the cases of a — 4, 5.2 and 8.5 
AU. To have a better insight on the difference in M cr i t among 
these cases, we analyse profiles for Z = 0.7 and Mp = 0.29M e 
for all the different cases shown in Fig.9. The results of changing 
a and the corresponding P out and 7' out (as dictated by Eqs.21, 22) 
are shown in Fig. 10 8 As expected from the behavior of M CI j t at 
Z = 0.7 in Fig.9, the cases of a < 3.2 AU follow the “dry” adia¬ 
bat, whereas the cases of a > 3.5 AU correspond to the “moist” 
adiabat. Since there are several factors that could be playing a 
role in shaping these profiles, we study the effect of changing 
just one boundary condition and fixing all the others with the 
standard values (i.e, those given in Tablel) in Figs.11 a,b. 


8 If migration were taken into account, this change in the boundary 
conditions would be the way to compute the internal structure of the 
planets in the case that immediate adaptation to the disc temperature 
and pressure values were assumed. 


Fig. 11a shows the effect of changing the outer pressure, P out . 
That is, the P out used to solve the internal structure equations is 
the one given by Eq.22, but all the other boundary conditions 
correspond to the standard case. From this figure we infer that 
the critical core mass is reduced when P out increases, since the 
pressure in the envelope is always higher for higher / J out . This 
fact can be understood as follows: from Eqs. 10, 11 and 13; it 
can be inferred that (|£)-/- < 0. Hence, at the outer boundary of 
the planet, the entropy decreases for larger P out . For these pro¬ 
files of Z=0.7, the envelopes are fully convective, so the entropy 
of the envelope is fixed by P out , 7’out. On the other hand, the tem¬ 
perature at the surface of the core is the same for all the cases 
because they correspond to the same planetary mass 9 . There¬ 
fore, the pressure at the surface of the core also has to be larger 
for larger P out (see Wuchterl 1993; Ikoma et al. 2001, for a dis¬ 
cussion on the sensitivity of the critical core mass with the outer 
boundary conditions for fully convective envelopes). 

In Fig.lib, the only parameter we change is 7'„ ut . The first 
remarkable observation is that if P out is larger than the temper¬ 
ature given by the Clausius-Clapeyron relation at Pq, no con¬ 
densation of water occurs. In our particular disc toy model this 
happens for a < 3.5 AU. Then, the change in 7’out is what gives 
raise to the two types of profiles observed in the general case 
(Fig. 10, where all boundary conditions are modified consistently 
with Eqs.21, 22). The general trend of changing P out is that M cr j t 
diminishes when P out diminishes. This can be explained draw¬ 
ing similar arguments than with Fig.1 la. Due to Eq.13, we know 
that (|f )/> > 0 , and hence, at the outer boundary of the planet 
the profiles with lowest 7’out are the ones with lowest entropy. 
Since the temperature at the surface of the core is the same for 
all the profiles, the profiles with lowest entropy will have the 
largest central pressure (because (ff )r is always negative). Con¬ 
sequently, regarding 7’out, the further the planet is from the star, 
the smaller will be M cr i t . 

We have found, therefore, that the change of P out and 7’out 
provoke opposite effects in the change of the critical core mass 
as we move further out in the disc. This is what makes Fig.9 and 
Fig. 10 not so easy to interpret. Fooking back at the profiles of 
Fig. 10, we see that the case of a =3.5 AU is the only one in which 
the profile is moist due to the increase in P out (when compar¬ 
ing to the reference case), which place the pressure-temperature 
curve in the left side of the Clausius-Clapeyron curve. For the 
cases of a =2, 3 and 3.2 AU, the increase in P out that actually 
occurs with respect to Fig.lib is not enough to switch the pro¬ 
files from dry to moist. Hence, in these cases, what makes these 
profiles follow the “dry” adiabat (i.e, water does not condense) 
is the increase in 7’out with respect to the reference case. 

Concerning the cases of a = 4, 5.2 and 8.5 AU, since 
(7’out,7’out) reside always in the region where H 2 O is in solid 
phase (or in the solid-vapor transition), the profiles follow the 
moist adiabat. Due to the opposite effects on M cl -j t that arise when 
Pout and Pout decreases as the planet is further from the star, the 
profiles in these cases are very similar. 

We can conclude then that the decrease in M cl i t with increas¬ 
ing metallicity due to the condensation of water in the envelope 
of protoplanets is robust with respect to changing boundary con¬ 
ditions. The exact location in the disc where H 2 O starts to con¬ 
dense will depend on the disc model and stellar spectral type, 
but for distances larger than this, the value of M cr j t will be prac- 


9 Property that can be inferred for stellar structure making use of the 
Virial Theorem, see, for instance, Kippenhahn et al. (2013). For bodies 
with a core this is not strictly true, but the changes in central temperature 
are negligible. 
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Fig. 9: Critical core mass as a function of metallicity for the 
specified semimajor axis and corresponding boundary condi¬ 
tions given by Eqs. 21 and 22. Dashed lines correspond to the 
reference cases (boundary conditions at a — 5.2 AU) plotted in 
Fig. 1: the yellow curve is the “restricted” case (i.e., without con¬ 
densation of water), and the blue one the “non-restricted” case, 
(where water condenses). For all the new cases (solid lines) wa¬ 
ter is allowed to condense. The only case where water does not 
condense for any value of Z is for a - 2 AU. 



T[K] 

Fig. 10: Pressure-temperature profiles for different position in 
the disc for Z-0.1 and Mp = 0.29M@. T out and P out correspond 
to those given by Eqs.21, 22 for the semimajor axis specified in 
the labels. 


tically independent of the location of the planet in the disc, as 
shows Fig.9 for a > 4 AU. 

6. Discussion 

6.1. Validity of the uniform metallicity assumption 

Although the disruption of planetesimals is more likely to oc¬ 
cur in the densest regions close to the surface of the core, the 
final fate of the “metals” is not a priori clear. The heavy ele¬ 
ments could sink slowly to the core, or get well mixed in the 
H-He envelope, depending upon the strength of convection and 
the solubility of the material in the H-He envelope. In Sect. 5.2 
we assumed that the metallicity was uniform throughout the en- 



T [K] 



T[K] 

Fig. 11: Pressure-temperature profiles for different position in 
the disc for Z-0.1 and Mp = 0.29M e . a) P out given by Eq.22 for 
the labeled semimajor axis, but 7' ou , = 150 K (reference case, see 
Table 1) for all the profiles . b) T out given by Eq.21 for the labeled 
semimajor axis, but P out corresponds to the reference case for all 
the profiles. 


velope, and this may not be always the case. Nevertheless, if the 
envelope is fully convective, this assumption is a reasonable one. 
This was the case for all envelopes with Z > 0.45. However, it 
could also be that the envelopes were fully convective because 
a uniform Z was assumed. It is important to know if this is the 
case or not, because if the enrichment does not reach the lowest 
temperature region of the envelope, no water will be present for 
temperatures lower than its condensation temperature, and then 
the condensation of water will not take place (at least for the 
boundary conditions considered in this work). 

To evaluate how well mixed the envelope will be, we per¬ 
formed the following test. We assumed that the enrichment of 
material occurs at temperatures higher than the condensation 
temperature of water. Hence, the outermost part of the atmo¬ 
sphere (T < T c end) is assumed not to be enriched and have a 
solar composition, whereas the innermost part [T > 7' coniJ ) is as¬ 
sumed to have a uniform Z (higher than solar). 10 If with this two 
layered model atmosphere the envelopes are still fully convec- 


10 This scenario could take place if the water droplets formed where 
T < T cond precipitated, which is the reason why we chose r cond as the 
boundary between the polluted and non polluted layers. 
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tive, then a global mixing of material will occur and therefore 
the envelopes will acquire, eventually, a uniform composition. 

Since this initial set-up implies a considerable jump in mean 
molecular weight at T - 7' coniJ (where the boundary between 
the solar composition and the high-Z region occurs), the Ledoux 
criterion for stability against convection should be used (Kippen- 
hahn et al. 2013, chap. 6 ), instead of the Schwarschild criterion. 
Therefore, for our calculation the radiative zones are those for 
which: 

V ra d < V ad + (23) 

where V„ = f , = (S) p j and 6 = ~(S For an 

ideal gas = 6 — 1 . 

Due to the finite jump of p at T = 7' con( |, V ;l diverges at this 
temperature. Hence, to compute the composition gradient, we 
assumed the same simple prescription used in stellar evolution 
models to treat the problem of overshooting (Alibert et al. 1999), 
i.e, we assumed that the change of p is effective within a distance 
0.27/ of the radius at which T = r cond (H is the scale height of 
the atmosphere at T — 7 mniJ ), and we take V /( constant in this 
interval. 

For the opacity of the solar composition layer, we use, for the 
sake of consistency, the same opacities that we used before, but 
with a Z = Z Q for the Ferguson gas opacities. It is important to 
remark that for the range of temperatures of this (predominatly) 
H-He region, the Ferguson opacities are extrapolated, and there¬ 
fore, its values are probably not very accurate. We will discuss 
the effect of varying this opacity in the next section. 

The results of the test show that the critical envelopes of 
the standard, non-restricted case remain fully convective when a 
layer of solar composition is considered for T < 7’ cond , provided 
that Z > 0.6. This means that for these high enriched cases, the 
assumption of uniform Z was a good one, and hence that con¬ 
densation of water is likely to take place, at least for the opacities 
considered here. 

Finally, whatever the fate of the water droplets (this depend¬ 
ing on the dynamics of the atmosphere and the miscibility of the 
forming droplets), if the envelope remains convective up to the 
layer where H 2 O can condense, it is likely that the continuous 
saturated state of H 2 O can sustain a moist-adiabatic pressure- 
temperature profile . 11 This will make it possible for the critical 
core mass to reduce considerably, as shows Fig.l. Therefore, for 
high enriched envelopes (which should not be so rare during the 
formation scenario, as shown by Fortney et al. 2013), the effect 
of condensation of water can really play a role. 

6.2. Role of microphysics 
6.2.1. Ideal gas assumption 

Since in this work an ideal gas equation of state was assumed for 
the structure calculation, it is relevant to have an estimate of the 
error introduced by neglecting the degeneracy of the gas at high 
pressure. For this reason, a comparative test for H-He envelopes 
implementing CEA and SCVH EOS was performed. 

The results of the planetary mass as a function of the core 
mass (obtained by solving the internal structure equations) are 

11 Actually, the condition of a fully convective envelope to have perfect 
mixing is overly severe. In reality, vertical mixing tends to take place 
in regions well above the convective-radiative boundary. On Earth, 
namely, the convective region extends for 10 km over the planet sur¬ 
face, but vertical, efficient mixing of material occurs up to 80 km height 
(Goody 1995). 



Mcore 1^©1 

Fig. 12: Planetary mass as a function of the core mass given 
by solving the internal structure equations for pure H-He en¬ 
velopes. The red curve uses CEA for the EOS and the blue, 
SCVH. Boundary conditions are those from Table 1. AM cr it ~ 
1.4 M e 


shown in Fig. 12. From this figure, it can be appreciated that dif¬ 
ferences in planetary mass for a given core mass start to be no¬ 
ticeable when reaching the critical core mass. The critical core 
masses are 14.18 M e and 15.58 M@, for CEA and SCVH cases, 
respectively; which implies a relative error of ~ 10% in the com¬ 
putation of the critical core mass. Even performing the same test 
with other boundary conditions or opacities, we find that the er¬ 
ror is never larger than ~ 12% 12 

To have a better insight of the effect of degeneracy, we show, 
in Fig. 13, a map of the difference in density obtained using CEA 
and SCVH EOS, as a function of (log P, log T). Overlapped are 
the profiles of the envelope obtained for Mp = M^(M rnr C EA) 
for CEA and SCVH cases (same conditions as in Fig. 12 were 
assumed). Although the difference between CEA and SCVH 
profiles is not easily distinguishable, it can be noticed that the 
CEA profile overestimates the central temperature and pressure, 
leading, as a consequence, to an also overestimated central den¬ 
sity. The difference in central density between these CEA and 
SCVH profiles is of 0.1 g/cm 3 . The higher density obtained us¬ 
ing CEA translates into a higher envelope mass, and hence, in 
a smaller critical core mass, as shows Fig. 12. We have shown 
that the change in critical core mass induced by considering the 
enrichment of the envelope (and in particular, the enrichment 
by condensable species) is orders of magnitude larger than the 
changes produced by including non-ideal effects. This justifies 
our assumption in assessing the effect of condensation on the 
determination of the critical mass. 


6.2.2. Remarks on the opacity 

Concerning the opacity, it should be kept in mind that the Fergu¬ 
son opacities for the gas are extrapolated in the coldest layers of 
the envelopes computed in this work. In principle this does not 

12 Note, however, that due to the non-linearity of the planet formation 
process, this difference in critical core mass can significantly influence 
the final characteristics of the planet. 
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Fig. 13: Map of the density difference between CEA and SCVH 
EOS as a function of temperature and pressure, for pure H-He 
composition. Superposed are shown the profiles of the envelope 
for the two different EOS for M v = 18.7 M 0 (planetary mass for 
which M core = M cr it when CEA is used for the EOS). All units 
in cgs. 



Fig. 14: red line: profile for two-layered atmosphere, but using 
Freedman opacities in the H-He region instead of the Ferguson 
opacities used in all other results of the paper. Blue line: Z=0.7 
(and uniform in the envelope) where water is allowed to con¬ 
densed. The resemblance of both profiles, compared to the “re¬ 
stricted” profile of Fig.5, suggests that for the red curve the cor¬ 
responding critical core mass will also be very low. Solid lines 
stand for convective regions and dotted for radiative ones. 


affect the profile of envelopes with high and uniform metallicity, 
since they are fully convective; but it could affect the structure 
for the test described in Sect.6.1. For instance, if instead of us¬ 
ing the extrapolated Ferguson opacities for the H-He region, we 
used Freedman et al. (2008) tables, the envelope would become 
radiative in the H-He region. This is because the Freedman opac¬ 
ity values are lower for the cold regions than the ones extrapo¬ 
lated from Ferguson, and hence, the adiabatic gradient becomes 
larger than the radiative one. In this case, then, the mixture of the 
material would probably not occur in the whole envelope. How¬ 
ever, because of the use of reduced opacities in the outermost 
part of the atmosphere, the critical core mass would still remain 
very low (Ikoma et al. 2000; Hubickyj et al. 2005; Hori & Ikoma 
2010), as suggested by Fig. 14. It is interesting to remark that 
even if condensation of water does not take place, the critical 
core masses shown in Fig. 1 will not be significantly affected. 

6.3. Chemical equilibrium assumption 

There are basically two processes that can preclude chemical 
equilibrium in planetary atmospheres: UV-photochemistry in 
close-in exoplanets, and vertical mixing in cold atmospheres of 
giant planets (Moses 2014, and references therein). It is well 
known that in Jupiter, the measured CO abundance at the top 
of the atmosphere is much higher than expected from chemical 
equilibrium calculations (Lodders 2010). The main reaction pro¬ 
ducing CO in planetary atmospheres is H 2 O+CH 4 3 IT +CO. 
The timescale of this reaction is indeed so long for the temper¬ 
ature and pressure at the top of Jupiter’s atmosphere, that CO is 
no longer destroyed by this reaction (Seager 2010). The amount 
of CO measured at the top of the atmosphere corresponds to that 
at T~1000 K, where chemical equilibrium still holds (Lodders 
2010). The CO produced at this level is carried outwards by 
vertical mixing. Hence, the CO is froze-in, or “quenched” at T 
~1000 K (the so-called quench level). 

The aforementioned reaction affects as well the amount of 
water at the top of the atmosphere, issue of major importance for 
our results, since the effect of water condensation on the critical 


mass depends on how much water is present in the cold layers of 
the envelope. 

The real amount of water will depend on several mechanisms 
operating in a planet formation scenario: the amount present in 
the incoming planetesimals, the depth at which they are com¬ 
pletely destroyed, the strength of convection and the extent of 
the convective cells, etc. To simplify this complex scenario, let 
us imagine the following picture. We ignore where planetes¬ 
imals deposit their mass. We assume large scale-convection. 
Hence, the volatiles provided by planetesimals (whatever their 
molecule/atomic state) will also be transported to the top of the 
core, where temperature is high enough to dissociate any water 
molecule. Then, due to convection, the bubbles reach the quench 
level and the abundance of water will remain fixed in that amount 
up to the top of the atmosphere. Therefore, in order to test in an 
approximate way the effect of disequilibrium, we assume in the 
test of this section that the amount of water is the one in equi¬ 
librium at a quenching temperature. In reality, the quench level 
depends on the species present in the atmosphere and their con¬ 
centrations. Given the similarity of the P-T profiles in our planets 
and in Jupiter around T= 1000 K, we take 7 qU ench= 1P00 K. 

In this section we focus on the abundance of water, which 
is the key factor in our results. For this, we consider a fictitious 
atmosphere that is made only of water and noble gases. Since 
no reaction occurs between water and noble gases, this scenario 
provides a good opportunity to study in an approximate way the 
effect of water condensation in a system where all the species 
remain quenched until the temperature is high enough for water 
to dissociate (i.e, for T<2000 K). We mimic, in this way, what 
would happen in a disequilibrium scenario. 

In Fig. 15 we present temperature-pressure profiles for the 
“nominal” composition case (the envelope composition assumed 
so far in this work, i.e, H, He, C and O) and the new “simplified” 
case, where no chemical reactions take place since the envelope 
is made of water and noble gases 13 . Both cases share the same 

13 The noble gases considered for the “simplified composition” are He 
and Ne and the mole fractions are chosen to have approximately the 
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Fig. 15: Comparison of pressure-temperature profile between the 
simplified envelope composition and the nominal composition 
for Z=0.7. In both cases water is allowed to condense. The same 
mole fraction of water at the top of the atmosphere is assumed 
for both profiles. Therefore, the departure from the moist adiabat 
occurs at the same (P,T) in both profiles. 

mole fraction of water at the top of the atmosphere. As was ex¬ 
plained in Sect.5.2, this guarantees that the departure from the 
moist-adiabatic profile occurs at approximately the same (P,T). 
Notice, however, that the profiles follow different tracks after 
the departure from the moist adiabat. This is due to the lack of 
chemical reactions in the simplified case. In the nominal case, 
chemical reactions take place at T> 500 K. This lowers the adi¬ 
abatic gradient and causes the bump of the blue curve in that 
range of temperature. Hence, computing the critical core mass 
with the simplified envelope composition provides us not only a 
way of simulating a quenched (or disequilibrium) scenario, but 
also allows us to disentangle the effect of reduction of critical 
core mass due to chemical reactions and water condensation. 

The effect of water condensation is determined by the point 
in the temperature-pressure diagram at which the vapor pressure 
stops being saturated (and this depends just on the mole fraction 
of vapor). Thus, concerning the effect of water condensation on 
the critical core mass, the two profiles shown in Fig. 15 are equiv¬ 
alent. We can therefore compute the critical core mass for the 
simplified case in two steps: 

1. We compute the critical core mass as a function of the mole 

fraction of water (xh 2 o) for the simplified composition (Fig. 

16). 

2. We compute the mole fraction of water as a function of Z for 

the nominal composition (Fig. 17, purple curve). 

Hence, by combining Figs. 16 and 17 (purple curve) we can 
express the critical core mass of the simplified composition as 
a function of an “equivalent” Z, for which “equivalent” means 
having the same mole fraction of water as the one on the top 
of the atmosphere of the nominal composition. For instance, for 
a mole fraction of water of 0.29 in the simplified case, M„ j t « 
0.5M® (Fig. 16). For the nominal composition, this mole fraction 
of water at the top of the atmosphere corresponds to a metallicity 
of Z=0.7. Thus, the simplified composition with a mole fraction 
of water of 0.29 has an“equivalent metallicity” of 0.7. 

same mean molecular weight at the top of the atmosphere than in the 
“nominal” composition case. 


Fig. 16: Critical core mass as a function of the mole fraction of 
water in the simplified case (purple) and nominal case (blue). 
Water condensation takes place in both cases. 



Fig. 17: Nominal composition case. Relation between the metal¬ 
licity and the mole fraction of water at the top of the atmosphere 
(blue) and at T= 1000 K (green). 


Nevertheless, as was explained before, the fact of having 
xh 2 o=0.29 at the top of the atmosphere in the nominal case for 
Z=0.7 is a result of assuming chemical equilibrium in that layer, 
which is probably not true. Thus, a more realistic amount of wa¬ 
ter at the top of the atmosphere should be the one at T= 1000 
K (namely, xh 2 o=0.22 for Z=0.7, see green curve of Fig. 17), 
where chemical equilibrium is a reliable assumption. We there¬ 
fore compute again the critical core mass for the simplified com¬ 
position, but assuming that the mole fraction of water is the 
one of the nominal composition at T=1000 K. The critical core 
mass as a function of the equivalent metallicity for the simplified 
compositional cases (with and without water condensation) are 
shown in Fig. 18. 14 

From Fig. 18 we can draw the following conclusions. In a 
chemical disequilibrium scenario (like the one we mimic with 
our simplified composition), considering or not the condensa- 


14 For simplification, fully convective envelopes were assumed. Since 
this is the case for the nominal composition for Z > 0.45, we just show 
results for this range of metallicities, so comparison with Fig. 1 is suit¬ 
able. 
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form, it does not take into account its condensation. However, 
we can at least test whether the abundance of NH 3 in the at¬ 
mosphere would be sufficient for its condensation to matter. For 
this, we consider an atmosphere in which we add nitrogen to the 
already included H, He, C an O. We keep the same solar C/O 
as used throughout, and assume also a solar ratio of N/O. Us¬ 
ing CEA, we find that in the coldest layers of the envelope, the 
mass fraction of NH 3 formed is just ~ 5% . Therefore, even if 
all the ammonia formed could condense, this would be a minor 
contribution compared to the amount of condensed water. Even 
more important, if condensation of other species took place, the 
effect of reduction of critical core mass would be even larger, 
since the effect of release of latent heat explained in Sect.5.2 is 
not unique to water, but qualitatively general to all vapor-solid 
(or vapor-liquid) phase transitions. 


Fig. 18: Critical core mass as a function of the equivalent metal- 
licity for the simplified composition cases. Purple solid line: case 
where water condenses and the mole fraction of water is the 
same as the one at the top of the atmosphere of the nominal 
composition case (i.e, the one given by the purple line of Fig. 17). 
Green-solid: water condenses but the mole fraction of water cor¬ 
responds to that at T= 1000 K of the nominal composition (green 
line of Fig. 17). Orange-dashed: same as green curve but in this 
case water is NOT allowed to condense. 


tion of water (green and orange curves of Fig. 18, respectively), 
brings a difference in critical core mass of still an order of magni¬ 
tude for large Z, the same that was found considering chemical 
equilibrium (Fig.l). Of course, the absolute values of M cnt are 
larger in the disequilibrium scenario, because the reduction of 
adiabatic gradient we had for T> 500 K (Fig. 5) due to chemical 
reactions does not take place in these water/noble gas envelopes. 
With this simplified composition we are excluding all chemical 
reactions (except dissociation of H 2 O). This is not justified for 
T> 1000 K, and also, in reality, the presence of dust grains may 
increase the rate of chemical reactions (Fegley & Lodders 1994) 
throughout the envelope. Therefore, the absolute A7 cn i values of 
Fig. 18 are upper limits when translating these into the nominal- 
composition out-of-chemical equilibrium case. 

Note finally that we assumed in this test that all water con¬ 
denses when the pressure and temperature are suited for it to 
occur. This is true in an equilibrium scenario. It is well known 
within the cloud formation community that in the presence of 
only one condensible species, the creation of the first conden¬ 
sates -nucleation- might require a high level of supersatura¬ 
tion (see, e.g., the reviews of Helling et al. 2008; Marley et al. 
2013). However, if the atmosphere contains preexisting nucleat¬ 
ing seeds, like dust grains, condensations can already occur for 
saturation ratios close to one. This heterogeneous nucleation is 
what explains cloud formation on Earth. In the atmospheres we 
are referring in this work, dust grains are much more abundant 
than on Earth due to the constant replenishment by the ablation 
of planetesimals. Therefore, water can effectively nucleate and 
condense in the atmospheres of these growing planets. 

6.4. Condensation of other species 

Apart from water, other compounds could condense in the 
pressure-temperature regimes considered in this work. In par¬ 
ticular, it is known that the atmosphere of Jupiter possesses NH 3 
clouds. Unfortunately, while CEA does handle NH 3 in gaseous 


6.5. Gas accretion timescale 


The growth of an enriched, critical protoplanet is not easy to 
infer. A simple estimation of the gas accretion timescale can be 
obtained by computing the Kelvin-Helmholtz timescale, which 
corresponds, approximately to (Ikoma et al. 2000 ): 


(r MM. 


tkh 


core iKZ env 




(24) 


Using Eq.7, this can be simplified to: 


M P1 


tkh 


Me, 


(25) 


Evaluating this expression when the structures become crit¬ 
ical, and using M core = 10 6 M, £l /yr (which is the approximated 
planetesimal accretion rate for the cases presented in this paper, 
as well as a standard value in the planet formation literature, 
see e.g., Bodenheimer & Pollack (1986); Pollack et al. (1996); 
Ikoma et al. (2000)), we find, for Z=0.7, that tkh ~ 7 x 10 4 yr 
when water condensation takes place, and tkh ~ 10 6 yr when 
water does not condense. In principle, these timescales would in¬ 
dicate that the condensation of water significantly speeds up the 
accretion of gas. However, two things have to be kept in mind: 


1. Accretion rate of gas depends on both the self-gravity of 
the envelope, and its ability to cool (and thereby, to con¬ 
tract). The latter depends on the microphysics of the enve¬ 
lope. Solving the full planet-growth problem for enriched 
envelopes would imply, in part, considering self-consistent 
grain growth calculation for the opacity (Ormel 2014; Mor- 
dasini 2014), because this affects the planet cooling. This, 
coupled to a code that solves the structure with a varying 
luminosity with time, will be implemented in a future work. 

2. Even if the gas accretion timescale is shorter for the case 
where water condensation occurs, gas accretion would lead 
to a decrease of the envelope metallicity (i.e. dilution), and 
therefore to the increase of the critical core mass. 


Indeed, if a planet reaches the critical mass at a given Z cnt , 
two effects will appear. First, as a result of solid (i.e., planetes¬ 
imal) accretion, the metal content of the planet (Mz) increases. 
Second, as a result of gas accretion, the mean Z of the planet 
will change. According to Ikoma et al. (2000), the gas accretion 
rate at the critical point is similar, but slightly larger, than the 
planetesimal accretion rate. At the time when the critical core 
mass is reached, the accretion of solids and gas will lead to a 
decrease of the mean Z in the envelope. As a result, the critical 
core mass, as well as the total amount of heavy elements for the 
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Fig. 19: Sketch of the possible phases an enriched protoplanet 
could follow after reaching the critical mass. Although not il¬ 
lustrated (for not overloading the image), planetesimal accretion 
is present in all the stages and the material assumed to remain 
in the envelope. Despite the competition between gas and plan¬ 
etesimal accretion, the immediate result would be dilution, since 
at this point the accretion of gas is larger than that of solids, as 
explained in the main text. As a consequence of dilution, M cr ; t in¬ 
creases. After this point, depending on the accretion rate of plan- 
etesimals (which gives a new Mz) versus the accretion rate of 
gas from the disc, the planet could become sub-critical or super¬ 
critical, as is discussed in Sect.6.5. 


critical structure Mz, c ut, increase. A key point, for the immediate 
future of the planet is therefore to compare the increase of Mz 
(resulting from the accretion of planetesimals) and the increase 
in the critical amount of heavy elements. 15 

Thus, two possible scenarios could take place after the crit¬ 
ical core mass is reached (see sketch on Fig. 19 for an illustra¬ 
tion). In the first scenario (case 1), the increase of critical mass is 
more rapid than the increase of Mz . In this case, the planet will 
become again subcritical, will not accrete anymore gas. Finally, 
as a result of planetesimal accretion, the metallicty may increase 
again. In the second scenario (case 2), the increase of critical 
mass is smaller than the increase of Mz . In this case, the planet 
will remain super-critical, and might become a giant planet. 

The occurrence of case 1 or 2 depends therefore on the rel¬ 
ative increase of M /cnl (Z) and Mz(Z). The steeper the slope of 
Mz, cnt (Z), the larger the increase of Mz, a it for a given reduction 
of Z. Thus, the higher the probability to be in case 1. For in¬ 
stance, as shown in Fig.20, when comparing the case where wa¬ 
ter condenses with the one where water does not condense, the 
former case has a steeper slope than the latter. Therefore, wa¬ 
ter condensation would lead probably to case 1, and hence to an 
increase of the heavy element content of the planet. 

In reality, which case takes place depends on the location of 
the forming planet, solid accretion rate, among other conditions. 


15 We are assuming in this discussion (as is illustrated in Fig. 19), that 
the mass of the core is fixed. Hence the heavy elements provided by 
planetesimals accretion remain in the envelope, making Z, and thereby 
M z , crit, grow. 



Fig. 20: Critical mass of metals (mass content of metals in the 
planet when M core = M cr jt) as a function of metallicity. Red line: 
restricted case. Blue line: non-restricted case. Boundary con¬ 
ditions given by Table 1. Note that despite the similarity with 
Fig.l, in this case the y-axis represents all the metal content of 
the planet at criticality, and not just that of the core. Depending 
on the slope of Mz, a it(Z), for the same planetesimal accretion 
rate, the planet could become sub or super-critical. The steeper 
the slope (water condensation case), the more likely is the struc¬ 
ture to become sub-critical, and thereby, to follow case 1 illus¬ 
trated in the Fig. 19. 


7. Conclusions 

The critical core mass for envelopes enriched in carbon and oxy¬ 
gen through the disruption of planetesimals was computed, cor¬ 
roborating the reduced values found by Hori & Ikoma (2011). It 
was found that large amounts of water are formed in the envelope 
of the enriched protoplanets, and that an important fraction of it 
can condense for boundary conditions suited for a >4 AU. Be¬ 
cause of the release of latent heat that takes place when H 2 O con¬ 
densation occurs, the energy budget of the envelope is dramat¬ 
ically changed. This translates into a drastic pattern shift in the 
pressure-temperature envelope profile, which implies a remark¬ 
able increase in pressure for a given temperature of the envelope. 
Therefore, compared to the case where H 2 O condensation is not 
included, the density of the envelope in the condensed case in¬ 
creases significantly, reducing the critical core to extremely low 
values, as shows Fig.l. 

The reduction of critical core mass due to water condensation 
depends on the amount of water present in the outer layers of 
the envelope, which is affected by the assumption of chemical 
equilibrium. However, we have shown that for a disequilibrium 
scenario, the critical core mass is still significantly reduced. 

This critical core mass is important in the context of the 
‘standard’ models of planet formation because it sets the time 
required to form a gas giant. Indeed, as shown in many papers, 
the gas accretion rate for a critical mass of ~ 10 M@ is so large 
that the formation time of a giant planet is for practical purposes 
very close to the time needed to reach the critical mass (see e.g. 
Bodenheimer & Pollack 1986; Pollack et al. 1996; Alibert et al. 
2005). 

However, this conclusion only holds for large values of the 
critical core, as the gas accretion rate at the critical mass is a 
rapidly increasing function of the planetary mass. As we have 
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shown here, the condensation of water reduces drastically the 
critical mass, and as a consequence, the gas accretion rate at the 
critical mass could be in this context very small compared to 
standard formation models. The work presented here shows that 
the standard picture of planet formation, including the concept 
of critical core mass, is overly simplified and should be revised 
in the future. 
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Fig. 4: Envelope composition as a function of temperature for Z = 0.2, 0.4, and 0.7. Mp=0.29 M® for all profiles (planetary mass 
when M core (Z=0.7) = M cnl (2=0.7)). Left column: non-restricted case. Right column: restricted case. In the non-resticted case it 
can be appreciated that water (solid lines) coexist in vapor phase (blue curve), solid phase (cyan curve) and even with liquid phase 
(violet curve) in the case of Z = 0.7. Except H 2 O in solid and liquid phase, the other 13 species shown here were readily considered 
in HI11. In the restricted case, water only exists in vapor phase. 
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